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VELOCITY AND TEMPERATURE DISTRIBUTIONS OF COAL-SLAG LAYERS 
ON MAGNETOHYDRODYNAMIC GENERATOR WALLS 
by Carlson C. P. Pian and J. Marlin Smith 
Lewis Research Center 

SUMMARY 

Approximate analytical expressions are derived for the velocity and temperature 
distributions in steady-state coal slag deposits flowing over magnetohydrodynamic gen- 
erator walls. Effects of slag condensation and Joule heating are included in the anal- 
ysis. The transport conditions and the slag temperature at the slag-gas interface are 
taken to be known parameters in the formulation. They are assumed to have been prede- 
termined either experimentally or from the slag properties and the gas dynamic calcula- 
tions of the free-stream flow. The analysis assumes a power-law velocity profile for 
the slag and accounts for the coupling between the energy and momentum conservation 
equations. Comparisons are made with the more exact numerical solutions to verify 
the accuracy of the results. 


INTRODUCTION 

The behavior of slag layers in magnetohydrodynamic (MHD) generator ducts has re- 
cently gained increased attention after the University of Tennessee Space Institutes f ex- 
periments showed the possibility of power generation from MHD generators coated with 
coal slag (ref. 1). Problems involved with analyzing these slag deposits are not 
straightforward because of the vast variation in physical properties and experimental 
conditions. Theoretical analysis, even with idealized assumptions, is complicated by 
the coupling of the slag momentum and energy equations through the strong dependence 
of the slag viscosity on the temperature. Hence, most investigators have relied on nu- 
merical procedures (refs. 2 and 3). Although analytical solutions provide more physical 
perceptivity, exact solutions for the velocity and temperature distributions within a slag 
layer have been obtained previously only by excluding slag condensation and Joule heat- 
ing effects (ref. 4). 

In the present report, approximate analytical expressions are derived for the veloc- 



ity and temperature variations in steady-state slag layers. Effects of mass condensation 
and Joule heating are included in the analysis. The transport conditions and the slag 
temperature at the slag-gas interface are taken to be known parameters in the formula- 
tion. They are either experimentally measurable or determinable from the slag proper- 
ties and the gas dynamic calculations of the free-stream flow. The analysis, based on 
the model of reference 5, assumes a power-law slag velocity profile and accounts for 
the coupling between the energy and momentum conservation equations. Comparisons 
are also made with the more exact numerical solutions of reference 4 to verify the ac- 
curacy of the present analysis. 


GOVERNING EQUATIONS 

Consider the steady slag flow over a smooth horizontal surface of infinite extent in 
the streamwise x-direction, as shown in figure 1. The assumption is made that the slag 
surface is acted upon by gas dynamic heating and shear stress at the slag surface t. 
from the concurrent gas flow. In addition, there may be mass transfer to the slag layer, 
at a rate m., from condensation of slag particles. The wall at y = 0 is maintained at 
temperature T . The value of T w is dictated by the design of the wall segment and its 
coolant passages. The slag surface at y = 6 is at temperature T^ and moves with 
velocity u^. The value of T. or 5 is defined by the surface viscosity and the slag 
properties (ref. 6) or is experimentally measurable (ref. 7). The interfacial shear and 
the amount of mass transport to the slag layer are determinable from a Reynolds anal- 
ogy treatment for a flat plate (refs. 2 and 6). From the same flat-plate, turbulent- 
boundary-layer analysis, the heat flux q is given by H(T o0 - T^). Typical values of 
pertinent slag properties including the heat transfer coefficient H are presented in ta- 
ble I. Hence, we consider T^,, T., or 6, and T w as given parameters in this 

formulation. 

The once-integrated slag layer momentum equation, neglecting inertia forces at low 
Reynolds number, is 


M — = — (y - 6) + r. + j B(y - 6) (1) 

ay ax 1 y 

4 2 

Taking j < 1x10 A/m and the characteristic values given in table I, the terms in 
equation (1) are found to be of the following magnitudes: 
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— (y - 6) ~ P xlO ^ pascals 
3x a 

_3 

r. ~ P xlO pascals 
j y B(y - 6) ~ 40 pascals 

where P a is the gas pressure and varies from about 4x10® pascals at the generator en- 
trance to 0. 5x10® pascals at the exit (ref. 6). Thus, in a typical base-loaded MHD 
generator away from the duct exit the inertial forces, the pressure gradient dP/dx, and 
the j x B body force within the slag layer are all much smaller in magnitude than the 
shear force 3/3y(p 3u/3y) and can be neglected in the momentum balance. Thermal 
gradients parallel to the surface are assumed to be much smaller than the normal grad- 
ients and are neglected in the energy equation. The governing equations of continuity, 
momentum, and energy, respectively, reduce to 


Jh + 

ix = o 

(2) 

dx 

3y 


du _ 

T i 

(3) 

3y 

fi(T) 



and 


k 



3T i 2 
= pcv — ~ J- 

3y c r 


(4) 


where the thermal conductivity k is taken to be independent of temperature and compo- 
sition. The viscosity coefficient p, being a strong function of temperature, is written 
as 



where the exponent /3 has a value between 17 and 20 for typical slags in the tempera- 
ture range under study (ref. 6). The mass balance of the slag layer is given by 
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resulting from the integration of equation (2) between y = 0 and y = 6. 
A slag velocity profile of the form 



( 7 ) 


is assumed, where the exponent a is treated as an unknown, as in reference 5. This 
exponent is evaluated only after the expression for the temperature profile is obtained; 
thus, ot serves as a coupling between the energy and momentum equations. From equa- 
tions (2), (6), and (7), the normal component of the velocity in the slag layer is 


v = 



( 8 ) 


where m^ > 0 implies slag deposition. 


RESULTS AND DISCUSSION 


Non- Current- Carrying Slag Layers 


The temperature distribution in a non-current-carrying slag layer (i. e. , j = 0) can 
now be obtained. Substituting equation (8) into (4) gives 


d 2 T _ _ dT 

dy 2 k W dy 


The solution for the thermal gradient is 


dT _ q 
dy k 


exp 


-cm. /g +2 
k(a + 2)\ 6 a+l 


( 10 ) 


and q is the heat flux from the external flow to the slag layer at the interface. Expand- 
ing equation (10) in series form and integrating term by term yields the temperature 



distribution through the layer as 


k ( T _ T 

q 


w' 


= exp 


cnijS 

k(a i + 2) 


cm. 


y - 


T a+ 3 


k(a? + 2)(a! + 3) g(j+l 



(ID 


For typical values of m^, c, k, and 6 (see table I), we can accurately approximate the 
temperature distribution by 


T = T +SLZ 


w 


exp 


cnou6 


k(a + 2) 


Slag velocity can now be obtained from equations (3), (5), and (12) as 


,.+ r 

IX. T§ / 


T +M 


w 


exp 


cm. 5 


k(a + 2) 


drj 


( 12 ) 


(13) 


u = D( <T w +^exp 


cn^S 


k k( a +2) 


1/3+1 


_ ^3+1 

w 1 


(14) 


where 


D a 


T ; k 


M t Tf(|3 + l)q 


exp 


-cm. 6 

k(a + 2) 


(15) 


Equation (14) can be observed to satisfy the no-slip condition at the wall by setting 
y = 0. Another expression for D incorporating the interface slag velocity can be ob- 
tained by evaluating equation (14) at the interface where y = 6, giving 


Hence, 



exp 


cm.6 
k(« + 2) 


1/3+1 



(16) 
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U ss u 



+ 91 


exp 


cmjfi 



k (a + 2) 


exp 


cmjS 


)fi+l 


Ha + 2) 


rj^+1 

W 


J 


-1 


(17) 


If we evaluate the derivative of u at the slag layer surface from equation (17) and equate 
this result to that obtained from equation (7), we find with use of equation (12) that 


(fl + D(T t - T w )Tf 

T0 + 1 _ T0 + 1 

i w 


(18) 


The slag layer thickness 5 is determined from equation (12) and evaluated at y = 6, 
expanding the exponential for small arguments, as 


6 = - (T. - T ) 
' 1 w 

<1 


1 - 


cm. 


q(a + 2) 


< T i • V 


(19) 


For nr = 0 this result reduces to that obtained in reference 6. 
From equations (12) and (14), the interface slag velocity is 


u i = 


T i k 


PjTfq^ + 1) 


exp 


-cm^5 


Ha + 2) 


(if +1 - Tg, +1 ) 


( 20 ) 


Substituting values of u^, 6, and a from equations (18) to (20) into equation (7) will give 
the velocity profile of non-current- carrying slag layers as a function of slag properties. 

The velocity profiles for various values of a. are presented in figure 2. One sees 
that most of the streamwise slag transport occurs near the slag surface when a » 1. 
This corresponds to the condition T? + '*‘ » T^ + '*\ In this limit, equation (18) reduces 
to 

T. - T 

a “ 03 + 1) — » 1 (21) 

T i 

As previously noted, (3 has a value between 17 and 20. In the limit T w — T., equa- 
tion (18) reduces to 


6 



a “ 1 + /3 


( 22 ) 


T i' 


w 


T i 


1 


and the velocity profile becomes linear. Figure 3 presents ck as a function of T w for 
various values of T. . In figure 4 the analytical result obtained in this report is com- 
pared to the more rigorous numerical result of Crawford, et al. (ref. 4). The agree- 
ment is excellent. For the case presented in reference 4, the value of gives £ =* 13 
for equation (5). 


Current-Carrying Slag Layers 

The temperature variation within a current-carrying slag layer flowing over an in- 
finitely segmented electrode wall is dictated by 

!i?T + f!Vxf +1 ® + i! = o (23 

dy 2 k \6/ dy ok 

In the present analysis, an electrical conductivity of the form cr(T) = AT y has been 
assumed (ref. 6). Furthermore, as part of the model, ct(T) is replaced by a(y). The 
relation between T and y needed for the conductivity calculation is provided by equa- 
tion (12) (i. e. , the temperature distribution when j = 0). Thus, the conductivity model 
applies only to cases where the contribution to the temperature distribution by Joule 
heating is much smaller in magnitude than 


T + 
w 




exp 


cm. 6 
x 

k(a + 2) I 


2 

In general, this implies j must be less than 10 000 A/m if T s 1000 K and j must 

2 W 

be less than 5000 A/m if T w < 1000 K. Using the heat flux at the gas-slag interface as 
the boundary condition, the first integral of equation (23) is 
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— = expj 

dy 


-cm.y' 


a +2 


lk(a + 2)6' 


:0t +l| 


cm. 6 

exp — 

k((x + 2) 

L J 


{2 


kA 


r 4-=£Li 

J [k(« + 2)5 £¥+1 . 


T +S2Z 


w 


exp 


cm. 6 


ivr 


k(a + 2) 


dq 


(24) 


Expanding those exponential terms that are functions of y for small arguments and in- 
tegrating the result term by term gives 


T(y) - T W +2X 


exp 


cm.6 


k(a + 2) 


.2 

J exp 


-cm. 5 


k(o + 2) 


Aq(l - y) 


T w +a8 exp 


cmjS 


k(o: + 2) 


il-y 


j 2 k exp 


-2cm i 5 


k(cy + 2) 


T + 3I 


Aq (1 - y)(2 - y) 


w 


k 


exp 


cm.6 


k(a + 2) 


)2-y 


j 2 kT 2_r exp 


-2cm. 6 


k(cK + 2) 


(25) 


Aq^(l - y)(2 - y) 


The remaining terms in this series expression for the temperature variation are at 
least a factor cm^6/k(o? + 2) smaller in magnitude. Their contribution is negligible 
under typical base-loaded conditions. Shown in figure 5 are the temperature distribu- 
tions at various current densities for two different T w ’s and im = 0. The following 
two relations have been used to approximate the numbers used for electrical conductivity 
in reference 4: 


a = 6. 892X10' 42 T 13 (ohm-m) -1 , T > 1341 K 


and 


0 = 1. 296X10" 10 T 3 (ohm-m) -1 , T < 1341 K 


8 



Direct quantitative comparison with the case presented in Crawford, et al. (ref. 4) is 
not possible since the present electrical conductivity model is not applicable to such a 
low value of T w and high value of j. 

The slag layer thickness can now be estimated by evaluating equation (25) at y = 6. 
The exact solution for 6 from the ensuing expression is extremely difficult, if not im- 
possible, to obtain analytically. However, since the magnitudes of the terms involving 
Joule heating effects have already been assumed to be small, the method of successive 
approximations may be used. As a first approximation, the small Joule heating terms 
are neglected, resulting in expressions (18) and (19) for a and 6. A better approxi- 
mation for the slag layer thickness can then be obtained by substituting the first approxi- 
mation solutions for a and 5 into the terms arising from Joule heating effects. The 
resulting expression is solved for 5. The result can be put in the form 



The slag velocity from equations (3), (5), and (25) is 



9 



Employing a binomial expansion in the form (1 + z) n , truncated to two terms, gives 


u(y) = — M 


M 


l 1 


+ 1) A q 3 X 3 (l - y)(j 3 + 1) 


T + &lf +1 - tP +1 

w k ) w _ 


-i 2 * 2 Tw(Tw * a ~) 1 ' Y jV^ 


w 


Aq 3 X 3 (l - y) 


Aq 3 X 3 (l - y)(2 - y. 




. 2 , 2 


.3.3 


Aq X (1 - y)(2 - y)(/3 - y + 2) 


+ 2) L' k / 


(29) 


wherein X s exp |];m.6/k(cy + 2)]. The velocity at the interface maybe found bv 
evaluating equation (29) at y = 6, noting that at the interface « [T + (qXS/k:]^. 
Thus, the interface velocity may be written as 


u i “ u i0^ + n) 

where u iQ is the value of u. when j = 0, given by equation (20), while 


(30) 


j 2 k exp 


il S 


t-2cm.61 

(a + 2)J 


Aq (1 - y) 


/3(y 2 - y/3 - 4y + /3 + SJT. 1 y + (/3 + l)/3T 2 " y+/3 Tr 1_ ^ 


(2 - y)(/3 - y + 2) 


- (0 + DT w T- y + 03 + l)Ti -1 Tj" y (2 - y) -1 

Similarly, comparing the gradient of the velocity at y = 5 evaluated from equa- 
tion (29) with that evaluated from equation (7) results in 

a ■“ a q (1 + <£) 

where a Q is the value of a when j = 0 given by equation (18) and 


( 31 ) 


( 32 ) 
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j 2 k exp( 


$ a. 


-2cm^TA 

/3q ) 


Aq 2 (l - y) 


(/3T. + T w )T: y - T 2_y T i : 1 (2 - y) 


-1 


J3(y 2 - y/3 - 5y + 2/3 + 5)1^ + (0 + l^T^^T: 1 ^ 
(2 - y)(j3 - 7 + 2) 

. i 1 - -K)T f~ V - (2 - V)T W T 1 1 ^ - T^ y j (3; 

(2 - y)(T. - T w ) J 

Typical behavior of $ and 5 as a function of current density, at various T w ’s, is 
shown in figures 6 and 7. As the amount of Ohmic heating is increased, the trend is 
toward thinner slag layers with more linear velocity profiles. 


CONCLUDING REMARKS 

Approximate analytical expressions are derived for the velocity and temperature 
variations in steady- state slag deposits flowing over the magnetohydrodynamic generator 
walls. Effects of both slag condensation and Joule heating are included in the analysis. 

In addition to providing trends and physical insight, the derived results enable one to 
easily determine the internal properties of coal-slag layers once the transport condi- 
tions and temperature at the slag-gas interface are known from preliminary calculations 
and/or experimental measurements. Because of the assumptions in the electrical con- 
ductivity model, the results for electric- current- carrying slag layers are not applicable 
to low -wall-temperature and high-current-density cases. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 19, 1976, 

506-25. 
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APPENDIX - SYMBOLS 


A a constant 

B magnetic field strength, T 

c specific heat at constant pressure of slag layer, J/(kg)(K) 

H heat transfer coefficient, qCT^ - T^ -1 , W/(m 2 )(K) 

p 

j,jy normal electric current density, A/m 
k thermal conductivity of slag, W/(m)(K) 

n^ mass transfer rate at slag surface, kg/(m 2 )(sec) 

P pressure, Pa 

P„ gas pressure, Pa 

q heat transfer rate to slag surface from free -stream gas, W/m 2 
T temperature, K 

temperature at slag surface, K 
T w wall temperature, K 

free-stream gas temperature, K 
u velocity parallel to wall, m/sec 

u. velocity o c . slag at y = 6, m/sec 

v velocity perpendicular to wall, m/sec 

x coordinate parallel to wall 

y coordinate perpendicular to wall 

a exponent in velocity profile 

/3 exponent in slag viscosity model 

y exponent in slag electrical conductivity model 
6 slag layer thickness, m 

2 

i u slag viscosity, N-sec/m 

viscosity at slag surface, N-sec/m 2 

O 

o slag density, kg/m 

a electrical conductivity of slag, (ohm-m)"* 

p 

r. shear stress at slag surface, N/m 
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TABLE I. - SOME TYPICAL SLAG PROPERTIES AND PARAMETERS 


IN BASE -LOADED GENERATOR DESIGNS 


Parameter 

Typical value 

Reference 

Thermal conductivity of 

0. 5 to 1. 73 W/(m)(K) 

3,6 

slag, k 

2. 5xl0 3 to 3. OxlO 3 kg/m 3 


Slag density, p 

2, 3,6 

Specific heat of slag, c 

1. 382X10 3 N-m/(kg)(K) 

2 

Slag layer thickness, 6 

0. 5xl0 -3 to 2. 5X10' 3 m 

2,6 

Slag layer Reynolds num- 

10" 1 

2 

ber 

1. 43x10 3 W/(m 2 )(K) 


Heat transfer coeffi- 

3,4 

cient, q/(T oo - T.) 



Rate of slag deposit, m. 

3. OxlO" 3 to 2. 8X10" 2 kg/(m 2 )(sec) 

6 

Shear stress at slag sur- 

P a xl0“ 3 pascals 

6 

face, Tj 



Interface velocity, u. 

~0. 4 m/sec 

2,4 

Gas- slag interface tem- 

1550 to 2000 K 

2,6 

perature, T. 



Wall temperature, T w 

500 to 2000 K 

2,4,6 

Free -stream gas tem- 

2700 to 2800 K 

2,4,7 

perature, T^ 



Free-stream gas pres- 

0. 5x10 3 to 4x10 3 pascals 

6 

sure, P a 



Magnetic induction, B 

4 to 5 T 


Pressure gradient, dp/dx 

-P a xl0 -1 pascals 

6 






05 



1000 1500 2000 2500 3000 

Wall temperature, T^, K 

Figure 3. - Velocity profile exponent as a function of wall 
temperature at various slag surface temperatures. Ex- 
ponent in slag viscosity model, (3, 20; normal electric 
current density, j, 0. 


Slag velocity, m/sec 



Distance from wall in slag layer, mm 


Figure 4. - Comparison with results of Crawford, etal. (ref. 4). 
Thermal conductivity of slag, k, 1.73 W/(m)(K); interface ve- 
locity, 0.36 m/sec; normal electric current density, j, 0; 
free- stream gas temperature, T m , 2700 K; wall temperature, 

T w , 500 K; temperature at slag surface, Tj, 2075 K; heat trans- 
fer coefficient, H, 1430 W/(m*)(Kh mass transfer rate at slag 
surface, mj, 0. 



Distance from wall in slag layer, mm 



Figure 5. - Temperature distribution in slag layer for various current densities at wall temperatures 
of 500 and 1000 K. Thermal conductivity of slag, k, 1.73 WMmKK); free-stream gas temperature, 
T®, 2700 K; heat transfer coefficient, H, 1430 WJ(nj9(Kh mass transfer rate at slaa surface, m*,, 0; 
electrical conductivity of slag, o, 6. 8921 47x10" 42 T 13 at T> 1341 K and 1.296xl0 10 Tat T< 1341 K. 


.035 


.030 

.025 

.020 

.015 

.010 

.005 


0 


Wall temperature, 



Normal electric current density, j, A/nr 

Figure 6. - Typical behavior of $ as function of 
current density. Exponent in velocity profile, 
a = oqII + $); thermal conductivity of slag, k, 
1.73 WMmKK); free-stream gas temperature, 

T mi 2700 K; heat transfer coefficient, H, 

1430 W/lm 2 )(K); mass transfer rate at slag sur- 
face, nij, 0; exponent in slag viscosity model, 

P, 20. 



Slag layer thickness, 6, mm 



Figure 7. - Typical behavior of slag layer thickness as 
function of normal electric current density. Thermal 
conductivity of slag, k, 1.73 W/(m)(K); free- stream 
gas temperature, T*. 2700 K ; heat transfer coefficient, 
H, 1430 WMrrrKK); mass transfer rate at slag surface, 
mj, 0; exponent in slag viscosity model, (3, 20; slag 
layer thickness at j = 0, 6 q, 3 mm. 
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